What happens to Q-baUs if Q is so large? 
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In the system of a gravitating Q-ball, there is a maximum charge Qmax inevitably, while in 
flat spacetime there is no upper bound on Q in typical models such as the Affleck-Dine model. 
Theoretically the charge Q is a free parameter, and phenomenologically it could increase by charge 
accumulation. We address a question of what happens to Q-balls if Q is close to Qmax- First, without 
specifying a model, we show analytically that inflation cannot take place in the core of a Q-ball, 
contrary to the claim of previous work. Next, for the Affleck-Dine model, we analyze perturbation 
of equilibrium solutions with Q ~ Qmax by numerical analysis of dynamical fleld equations. We flnd 
that the extremal solution with Q = Qmax and unstable solutions around it are "critical solutions" , 
which means the threshold of black-hole formation. 

PACS numbers: 04.25. dc, 04.40.-b, 11.27.-|-d, 95.35.-|-d 



I. INTRODUCTION 

In a pioneering work by Friedberg et al. in 1976 [l|, 
nontopologieal solitons were introduced in a model with 
a U(l)-symnietric complex scalar field coupled to a real 
scalar field. In contrast with topological defects, they 
are stabilized by a global U(l) charge, and their energy 
density is localized in a finite space region without gauge 
fields. In 1985 Coleman showed such solitons exist in a 
simpler model with an S0(2) [viz. U(l)] symmetric scalar 
field only, and called them Q-balls [2| . 

Q-balls have attracted much attention in particle cos- 
mology since Kusenko pointed out that they can exist 
in all supersymmetric extensions of the standard model 
0] . Specifically, Q-balls can be produced efficiently in the 
Affleck-Dine (AD) mechanism [J] and could be responsi- 
ble for baryon asymmetry [11 and dark matterQ . Q-balls 
can also infiuence the fate of neutron stars [31 • Since Q- 
balls are supposed to be microscopic objects, equilibrium 
solutions and their stability have been intensively stud- 
ied in flat spacetime Q. It was shown that catastrophe 
theory is a useful tool for stability analysis of Q-balls Q . 

If Q-balls are so large or so massive, on the other 
hand, their size becomes astronomical and their gravi- 
tational effects are remarkable [l3, [HI ■ Such gravitat- 
ing Q-balls, or Q-stars, are analogous to boson stars p^ . 
While Q-balls exist even in flat spacetime, boson stars are 
supported by gravity and nonexistent in flat spacetime. 
Multamaki and Vilja showed that the size of Q-balls is 
bounded above due to gravity ll| . Becerril et al. studied 
evolution of unstable solutions by numerical analysis of 
dynamical field equations [Tsj . 

In our previous work [l^ - [l6j we studied equilibrium so- 
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TABLE I: Comparison of maximum and minimum charge (or 
gravitational mass) for fiat Q-balls, gravitating Q-balls and 
boson stars. 



lutions and their stability for several models using catas- 
trophe theory to understand a comprehensive picture of 
flat Q-balls, gravitating Q-balls and boson stars. One of 
the main results is summarized in TABLE I. 

• In flat spacetime, as long as the absolute minimum 
of the potential V{4i) is located at (/> = 0, there is 
no upper bound on charge. If we take self-gravity 
into account, however, there is maximum charge 
Qmax (or maximum energy i^max), at which stabil- 
ity changes, regardless of models. 

• In flat spacetime, in some models such as the AD 
gauge-mediated model there is nonzero mini- 
mum charge Qmin, where Q-balls with Q < Qmin 
are nonexistent. If we take self-gravity into ac- 
count, however, there exist stable Q-balls with ar- 
bitrarily small mass and charge. 

The above properties of gravitating solitons hold for gen- 
eral models of Q-balls and boson stars as long as the lead- 
ing order term of the potential is a positive mass term in 
its Maclaurin series. 

Theoretically the charge Q is a free parameter, and 
phenomenologically it could increase by charge accumu- 
lation. Therefore, a natural question could arise: what 
happens to Q-balls if Q is close to Qmax? In this paper we 
shall investigate what happens if we give perturbations 
to equilibrium Q-balls with Q w Qmax- 
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In connection with this question, Matsuda [l7[ claimed 
that inflation occurs in the core of a Q-ball if Q is large 
enough. He assumed a kind of hybrid potential. In his 
scenario inflation takes place when the gravity-mediated 
term, 
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dominates, where m is the gravitino mass, K term a one- 
loop correction, and M the renormalization scale. If this 
scenario is true, it would have an important implication 
for inflationary models. Therefore, the present study is 
also important as a close examination of this scenario. 

This paper is organized as follows. In Sec. II we review 
previous results of gravitating Q-balls in the AD gravity- 
mediated model. In Sec. Ill we make analytic discussions 
on general properties of Q-ball equilibrium solutions. In 
Sec. IV we discuss whether Q-ball inflation can take 
place. In Sec. V we present dynamical field equations 
and our computing method. In Sec. VI, by analyzing the 
dynamical field equations, we investigate what happens 
if Q is so large in the AD gravity-mediated model. Sec. 
VII is devoted to concluding remarks. 



II. GRAVITATING Q-BALLS IN THE 
AFFLECK-DINE POTENTIAL 

To begin with, we review previous results of equilib- 
rium solutions of gravitating Q-balls in the AD gravity- 
mediated model (|l.ip [l5| . Consider an S0(2)-scalar field 
4> = {4>i,4>2) coupled to Einstein gravity. The action is 
given by 
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(2.1) 

where (/) = cj)\+ cj)\. We assume a spherically symmet- 
ric and static spacetime. 



and homogeneous phase rotation, 

cf){t,r) = (f){r) {cos utt, sin ujt). 
Then the field equations become 
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FIG. 1: Existing domain of equilibrium solutions for K = 
—0.01 in the Q-E space. We choose «: = 0, 0.01. 
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The boundary conditions are given by 

0'(O) = 0(cx3) = A'{0) ^ a'(0) ^ 0, 

^(0) = a(oo) = 1. (2.7) 

The charge and the Hamiltonian energy are defined by 
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where Ms is the Schwarzschild mass. To perform numer- 
ical calculations, we rescale the relevant quantities as 
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K = GM^, 

(2.9) 



Figure 1 shows existing domain of equilibrium solu- 
tions for K = —0.01 in the Q-E space. In the case of 
flat spacetime (k = 0) no upper bound on Q ov E ap- 
pears and all the solutions in this flgure are stable |18| . 
In the case of k = 0.01, on the other hand, there appears 
a cusp A, which corresponds to Qmax and iSmax- The 
lower branch represents stable solutions, while the upper 
branch unstable ones. 

In flat spacetime equilibrium solutions exist only if 
K < 0. If we take self-gravity into account, however. 
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FIG. 2: Existing domain of equilibrium solutions for K = Q 
in the Q-E space. We choose k = 0.01. 

equilibrium solutions also exist even \i K > Q [l^. The 
solutions for K = Q have been known as mini-boson 
stars [l2|. Figure 2 shows existing domain of these equi- 
librium solutions in the Q-E space. The result is sim- 
ilar to that for K < 0. There appears a cusp, which 
corresponds to Qmax a-nd -Emax: due to gravity. The 
lower branch represents stable solutions, while the up- 
per branch unstable ones. 

One may wonder what happens to equilibrium solu- 
tions if K is so large. In the case of topological defects, 
static solutions are nonexistent if the vacuum expecta- 
tion value of the Higgs field is close to the Planck mass, 
and the defects expand exponentially instead [l^. By 
analogy with this topological inflation, one may conjec- 
ture that Q-ball equilibrium solutions are nonexistent if 
K is larger than some critical value of order one. 

To examine this conjecture, for the case of K = 
—0.01, we analyze equilibrium solutions with k = 
0.1, 1, 10, 100, 1000 and 10000, too. We show the ex- 
isting domains of the equilibrium solutions in Fig. 3(a). 
Contrary to the above expectation, we find that equilib- 
rium solutions exist even if k 3> 1. Furthermore, Fig. 
3(a) indicates 

K-Emax = constant, (2.10) 

where -Emax denotes the maximum energy for each 
K. Because the Schwarzschild mass of a Q-ball is 
given by Ms = 2^ |ll|, the relation (|2.10p means 
GMs =constant. Therefore, as long as (|2.10p is satis- 
fied, we expect that static regular solutions can exist. 

To confirm this argument, we present one of the metric 
functions A(r) of the extremal solutions with E-a^ax hi Fig. 
3(b). All configurations of A(r) are virtually the same, 
which is consistent with (|2.10p . We therefore conclude 
that equilibrium solutions exist no matter how large k is. 
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FIG. 3: Dependence of equilibrium solutions on k. We 
choose K = -0.01 and k = 0.01, 0.1, 1, 10, 100, 1000 and 
10000. (a) Existing domain of equilibrium solutions, (b) The 
metric function A{f) of the extremal solutions for each k. 



III. GENERAL PROPERTIES OF 
EQUILIBRIUM SOLUTIONS 

Our interest is perturbation of equilibrium solutions 
with Q sa Qmax- Before analyzing their evolution, how- 
ever, it is important to understand general properties of 
equilibrium solutions. 

In the case of flat spacetime, the field equation is given 
by 

0" + -0' + a;20=^. (3.1) 
r acp 

This is equivalent to the field equation for a single static 
scalar field with a potential \{j = V — uj'^cfP /2. If one 
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FIG. 4: Interpretation of Q-ball solutions by analogy with a 
particle motion in Newtonian mechanics. 



regards the radius r as 'time' and the scalar amplitude 
(j){r) as 'the position of a particle', one can understand Q- 
ball solutions in words of Newtonian mechanics, as shown 
in Fig. m Equation describes a one-dimensional 

motion of a particle under the conserved force due to 
the 'potential' — and the 'timc'-dependcnt friction 
— 20'/r. If one chooses the 'initial position' 0(0) appro- 
priately, the static particle begins to roll down the po- 
tential slope, climbs up and approaches the origin over 
infinite time. Because of the energy loss by the friction 
term in Newtonian mechanics, one finds 



V >0 at r sa 0. 



(3.2) 



Dominance of the kinetic energy over the potential energy 
in the core is one of the important properties of Q-balls, 
as we shall discuss below. 

We can extend the above argument to gravitating 
Q-balls by redefining the effective potential as = 
In this case, because 'the potential 
of a particle' is 'time'-dependent, the above explanation 
in Newtonian mechanics cannot apply. However, because 
a{r) is an increasing function, the 'potential' —Vuj{(j>) de- 
creases as 'time'. This means that some of the 'mechan- 
ical energy' is lost by this variation of the 'potential' as 
well as by the friction terms. Therefore, we obtain the 
relation, 
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We have confirmed that the inequality (|3.3I) holds for all 
our numerical solutions. 

The condition p.3|) is important because from it we 
can draw a general conclusion that the Q-ball core is 
attractive as follows. Because the regularity condition 
at r = is given by (j2.7p . we can expand the metric 
functions as 



(r) = «(0)(l + $(r)), A{r)^l + f{r), 



(3.4) 



with the boundary conditions 

$(0) = $'(0) = /(O) - /'(O) - 0. (3.5) 

In the vicinity of r = we can expand the geodesic equa- 
tions and the Einstein equations up to first order of $ 
and /, which yields 
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Note that we have not introduced the weak-gravity ap- 
proximation. 

Equations p.6|) and (|3.7|) correspond to an equation of 
motion and Poisson equation, respectively, in Newtonian 
mechanics. They lead to $" > at r « 0. Furthermore, 
because of the boundary condition $' (0) = , one finds 
$' > at r « 0, that is, the core region is attractive. 



IV. IS Q-BALL INFLATION POSSIBLE? 

It was argued that inflation can take place in the core 
of a Q-ball if Q evolves with tim e by absorbing other Q- 
balls and becomes large enough . In the last section, 
however, we showed that Q-balls have attractive nature, 
which indicates that Q-ball inflation is improbable. Here 
we shall discuss this issue more explicitly. 

Inflation, or accelerated expansion of a local region, is 
defined by the following conditions. 

1. Some local region is well-approximated by the 
Friedmann-Lemaitre- Robertson- Walker metric. 



-dT^ + a(T)2{dx' + X^ide^ + sin^ dip^)}. (4.1) 



2. d^a/dr^ > in this region. 

3. The volume of this local region increases. 

If any of these conditions are violated, we cannot say that 
inflation takes place. 

The two expressions (|2.2p and (|4.ip look very different. 
In the case of topological inflation, however, the core of 
a topological defect is described by de Sitter spacetime. 
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which assures us of consistence between the two expres- 
sions. If Q-balls also trigger inflation, the core of the equi- 
librium solutions should be well-approximated by (|4.ip . 

Here, under the assumption that the first condition is 
satisfied for the equilibrium solutions, we discuss whether 
the second condition is satisfied. The Einstein equations 
for the doublet scalar field yield 
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Although the relation between the two coordinate sets, 
(t, r) and (r, x), is not given explicitly, we find adt = dr 
at the center {r ^ x ~ Then wc can rewrite Eq. ()4.3|) 

as 
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at r « 0. 



(4.4) 



Because oj^ (j)'^ / a{0)'^ > 2V from ([13]), we find d'^a/dr^ < 
in the core. Unless we give so large perturbation that 
\d(p/dT\'^ /V is more than double, d^a/dr^ remains nega- 
tive. 

If we chose initial values of <p{x) and d4>/dT{x) ar- 
bitrary, the right-hand side of (|4.3p might be positive. 
However, if the initial configuration is far from that of 
the equilibrium solution, such infiation cannot be called 
Q-ball inflation. 

We should also note that we have not specified Q and 
a potential type V{(j)). We therefore conclude that Q- 
ball infiation cannot take place by charge accumulation 
regardless of Q and a potential type. 



V. DYNAMICAL FIELD EQUATIONS AND 
COMPUTING METHOD 



Although we have shown that inflation cannot take 
place even if Q > Qmax, it is still unclear what happens 
in this case. To address this question, we have to solve 
dynamical field equations. Here, in preparation for this, 
we present dynamical field equations and our computing 
method. 

We consider a spherically symmetric and dynamical 
spacetime, 

ds^ = -a^{t, r)df + A^{t, r)dr^ + r^{d6^ + sin^ 6dLp^). 

(5.1) 

Introducing dimcnsionless auxiliary variables. 
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we write down the field equations derived from the action 
(1^ as 
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We have regularized the dynamical equations at the cen- 
ter, in the sense that all of them contain no diverging 
term like 1/r. 

As for initial conditions, we assume 



d(0,r) = {A =)0, 
and the perturbed field configuration. 



(5.7) 



0(0, r) - (Mr) + 5^ir), 0), 0(0, r) = (0, c^0o(r)), 

(5.8) 

where (f>o is an equilibrium solution and Scjj is a small 
perturbation. For definiteness we adopt 



64>{r) = S(j){0) exp 



L2 



(5.9) 



where L is a length parameter. We rescale it as L = rriL 
and set Z = 2 in Sec. VI. 

To obtain initial values of a (or A) and a, we integrate 
(|5.3|) and ()5.5p with respect to r. More precisely, to keep 
better precision, only at the initial time we introduce an- 
other auxiliary function as Y = A — 1 = ar^ and rewrite 
(lOl) as 
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Our dynamical field variables are a, a, 0, ^ and za. 
Among them the lapse function a is determined by in- 
tegration of (|5.5I) with respect to r at each time step. 
The rest of the dynamical variables are determined by 
integration of (|5.4p and (|5.6p together with 
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which is given by the definition (|5.2p . The Hamiltonian 
constraint (|5.3p is not solved except for the initial values, 
but it is used to check numerical accuracy of the above 
time- integration. 

To perform integration of the dynamical variables with 
respect to we discretize a space with a mesh with a 
equal size Ar, 



r, = (i-l)Ar, 1 = 1,..., AT, 



(5.12) 



and label a dynamical variable F{t,ri) as f/*"*. Here 
F(t,ri) represents a, a, 0, ^ and -ccr collectively. We 
choose N ~ 3001. A derivative with respect to is 
approximated as 



dF ^ 9F 

— 2t(^'-) = — — (^0 



9(r2) 
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For each F^*^ we integrate with respect to t by the second 
order Runge-Kutta method, or the so-called "predictor- 
corrector" method. 
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FIG. 5: Perturbation of a stable solution A. (a) Field con- 
figuration of the equilibrium solution and a perturbed initial 
configuration, which is given by (|5.8I) with S(j){0)/(f>o{0) —0.1 
or -O.f and L = 2. (b) Time- variation of (fi{t, f — 0). 



As for the boundary conditions of the center and the 
outer edge, we follow Hayley and Choptuik [g^] as follows. 
For (p, ^ and -ca at r = 0, we employ a "quadratic fit," 
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For 4 and zu at the outer boundary, we employ 
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To suppress numerical errors further, we apply numer- 
ical dissipation to </)„,^„ and -c:7„, following Hayley and 



(a) 




FIG. 6: Perturbation of an unstable solution B: 
S4>{0)/(l)o{0) = 0.01 and Z = 2. (a) Time- variation of 
(j>{t,f = 0). (b) Snapshots of the metric functions A and 
a at t = 0, 4, 5, 5.4. 



Choptuik [131 • After the next value is evaluated, 



we set 



^(t+Af) ^ ^(t+At) 
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4^(t)^+^(t)^). (5_i6) 



where e is an adjustable parameter in the range < e < 1, 
and we choose e = 0.5. 



VI. WHAT HAPPENS IF Q IS SO LARGE? 

In the analysis method devised in the last section, we 
shall discuss what happens to Q-balls with the poten- 
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<^0(O)/0o(O) = -0.01 and L = 2. (a) Time-variation of 50(O)/0o(O) = 0.01 and Z = 2. (a) Time-variation of 
<P{i,r = 0). (b) Snapshots of the metric functions A and f = q). (b) Snapshots of the metric functions A and 

a at i= 0, 40, 80, 120, 160. a at t = 0, 20, 30, 37. 



tial if Q is so large. Because ordinary Q-balls exist 
only UK < 0, we concentrate on the case oi K < 0. 
Specifically, we analyze dynamical field equations by giv- 
ing perturbations to equilibrium solutions A, B, and C 
in Fig. 3. 

First, we consider perturbation of a stable solution 
A. Figure 5(a) shows the field configuration of the equi- 
librium solution A and perturbed initial configurations, 
which is given by with 5(j){0)/(t)o{0) =0.1 or -0.1 and 
L = 2. (b) shows the time- variation of (j>{t, f = 0), which 
indicates that the field continues to vibrate around the 
equilibrium configuration. These results assure us that 
stable solutions indicated by energetics or by catastrophe 
theory are really stable. The mean values of <j>{t, 0) for 



the two solutions are slightly different from each other 
because Q is slightly changed by the perturbed field (50. 

Secondly, we consider perturbation of an unstable so- 
lution B. We give two types of perturbations. Figure 6 
shows the case of positive perturbation, (50(O)/(/)o(O) = 
0.01. (a) indicates that the field diverges in the center, 
(b) tells us that a approaches to zero and A diverges. In 
the coordinate system (|2.2p this behavior means a black 
hole is formed. Figure 7 shows the case of negative per- 
turbation, 6(p{0)/<j)o{0) = -0.01. We find that the Q-ball 
diffuses most of mass and charge but not all. It becomes 
a thick-wall Q-ball with much smaller charge. 

Thirdly, we consider perturbation of the extremal solu- 
tion C with Qmax- Again, we give two types of perturba- 
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FIG. 10: Perturbation of the extremal solution for k = 10000 
and K = -0.01. We choose 5<t){Q) / <f)o{Q) = 0.01 and Z = 2. 
The figure shows snapshots of the metric functions A and a 
at i=Q, 20, 30, 35.5. 



way as in Fig. 8. 

Thus our numerical analysis has provided eonfirmation 
of our analytic argument that inflation cannot take place 
in the core of a Q-ball. Furthermore, it indicates that 
the extremal solution and unstable solutions near it are 



critical solutions of black- hole formation [21|. In fact, this 
critical phenomenon was already found for mini-boson 
stars {K = in the model (|l.ll) by Hawley and Choptuik 
[20| . It is reasonable that Q-balls with A' < and with 
if = share the same property in the case that Q = 
Qmax, or gravitational effects are so large. 



FIG. 9: Perturbation of the extremal solution C: 
J0(O)/(;/)o(O) = -0.01 and L = 2. (a) Time- variation of 
cj>{t,f — 0) and ma,x[A{t,r)]. (b) Snapshots of the metric 
functions A at i = 0, 40, 80, 120, 160, 200. 

tions. Figure 8 shows the case of positive perturbation, 
(5(/)(0)/(/)o(0) = 0.01. Like the unstable solution B, this Q- 
ball collapses and becomes a black hole. Figure 9 shows 
the case of negative perturbation, (50(O)/<^o(O) = —0.01. 
The behavior in this situation is not analogous to that 
for B in Fig. 7. The Q-ball continues to oscillate without 
diffusing mass or charge. The above results are also seen 
even if we take other values of L and (50(O)/(/)o(O) as long 
as they are not so large. 

Finally, we ascertain that the dynamics is virtually un- 
changed even if we choose k = 10000. Figure 10 shows 
an example of the dynamical solutions, where except for 
K the parameters are the same as in Fig. 8. We see that 
the Q-ball collapses and becomes a black hole in the same 



VII. CONCLUDING REMARKS 

We have addressed a question of what happens to Q- 
balls if Q is close to Qmax- First, without specifying 
a model, we have shown analytically that the core of 
an equilibrium Q-ball has attractive nature and infla- 
tion cannot take place there. Next, for the Affleck-Dine 
model, we have analyzed perturbation of equilibrium so- 
lutions with Q w Qmax by numerical analysis of dynam- 
ical field equations. We have found that the extremal 
solution with Q = Qmax and unstable solutions around 
it are "critical solutions", which means the threshold of 
black-hole formation. 

Specifically, for initial data (|5.8|) with (|5.9p . a black- 
hole is formed ii 6(f> > 0. U Scj) < 0, a, black hole is 
not formed, and there are two types of evolutions. If 
the initial configuration is very close to the extremal Q- 
ball with Q = Qmax, the Q-ball continues to oscillate 
without diffusing mass or charge. In other cases, the Q- 
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ball diffuses most of mass and charge, and becomes a 
thick-wall Q-ball with much smaller charge. 

Further study is necessary to understand detailed be- 
havior of these critical solutions. It is also interesting 
to investigate how such behavior of critical solutions de- 
pends on models V((j)). 



numerical calculations were carried out on SX8 at YITP 
in Kyoto University. 
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